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Abstract: Identification of fractional order systems is considered from an algebraic point of 
view. It allows for a simultaneous estimation of model parameters and fractional (or integer) 
orders from input and output data. It is exact in that no approximations are required. Using 
Mikusihski's operational calculus, algebraic manipulations are performed on the operational 
representation of the system. The unknown parameters and (fractional) orders are calculated 
solely from convolutions of known signals. A generalized Voigt model describing a viscoelastic 
material is used to illustrate the approach. 



1. INTRODUCTION 

Fractional order models have gained increasing interest 
over the last years. Torvik and Bagley (1984) gave one 
of the first mathematical justifications for the use of such 
models for viscoelastic materials. However, fractional mod- 
els have been utilized for a wide spectrum of physical sys- 
tems including batteries, magnet-suspension systems, elec- 
trical circuits, as well as biological and chemical systems - 
to name but a few. Several examples, including fractional 
systems with distributed parameters, can be found in the 
books Oldham and Spanier (1974) and Podlubny (1999). 

In addition, researchers from different domains have given 
experimental evidence for the usefulness of fractional mod- 
els by identifying their parameters and fractional orders. 
Some of the recently suggested identification procedures 
use fractional state variable filters based on a known 
fractional order (Cois et al. (2001)), frequency response 
functions (Kim and Lee (2009)), or finite element methods 
for the approximation of fractional derivatives (Schmidt 
and Gaul (2002)). An overview on system identification 
for fractional models can be found in Malti et al. (2007). 
Nevertheless, most methods known to the authors rely on 
some kind of approximation. 

An algebraic approach was used in Fliess and Sira- Ramirez 
(2003) to identify parameters in ordinary differential equa- 
tions and in Rudolph and Woittennek (2008) for partial 
differential equations. The present contribution extends 
the method to linear fractional models, both, with lumped 
and distributed parameters. The approach allows for the 
identification of system parameters and fractional orders. 
It gives exact relations, in the sense that no approxi- 
mation is required. Furthermore, no assumptions towards 
commensurability of fractional orders and system stability 
have to be made. 

The basic idea of the method is to use an operational repre- 
sentation of a fractional model (in the sense of Mikusihski) , 



usually described by an equation relating a (known) input 
and output, and to eliminate at least the unknown non- 
integer powers of s, corresponding to fractional derivatives, 
by some algebraic manipulation. Here, the focus lies on ob- 
taining an operational equation, the expressions in which 
can easily be interpreted as functions of time. Unknown 
quantities are calculated from relations involving only con- 
volutions of known (input and output) signals. 

The paper is structured as follows. First, basic background 
regarding fractional derivatives and Mikusihski's opera- 
tional calculus is given in section 2. Then, it is discussed 
how to obtain an equation suitable for the identification of 
fractional orders and model parameters in linear fractional 
models. Homogeneous initial conditions are treated in sec- 
tion 3, inhomogeneous ones 1 in section 4. A generalized 
Voigt model for a viscoelastic material serves as an exam- 
ple. In section 5 the results are extended to a fractional 
distributed parameter system, the diffusion- wave equation 
of fractional order, followed by a brief conclusion. 



2. GENERAL DEFINITIONS AND NOTATION AL 
ASPECTS 

In this section the most important basic definitions used 
in this paper are revisited without being exhaustive on 
technical details. For a elaborate introduction into the 
mathematical background of fractional systems the reader 
is referred to Oldham and Spanier (1974) and Podlubny 
(1999) or for a well written brief summary to Gorenflo and 
Mainardi (1997). The calculus of Mikusihski is detailed in 
Mikusihski (1983). 



as 



1 Only classical definitions of initial conditions arc considered here 
Recent findings by N. Maamri and J.C. Trigcassou as well 
T.T. Hartley and C.F. Lorenzo (Trigeassou et al. (2012), Lorenzo 
and Hartley (2008)) show that these may not be sufficient, in general. 



2.1 Fractional derivatives 

One of the (or maybe the) most commonly used defini- 
tion^) for fractional derivatives is the one attributed to 
B. Riemann and J. Liouville. A fractional derivative of 
order 2 a £ R+ of a function / is defined as 

i d v r* f(c-) , , 

da, ctf^v, 



a+l- 



dcr, a^v 



(3) 



D? f{t) = { T{y v ~ a) dt Jo {t ~ ^ a+1 ^ 
j t f(t), 

(1) 

t > 0, where v — \ < a < ^ e N and T denotes 
the Gamma function. The inverse operation of fractional 
differentiation is fractional integration. For t > and 
a e R + it is defined as 

J?f{t) = ^ r -j\t-oY- 1 f{cj)do-. (2) 

An alternative definition has been introduced by M. Ca- 
puto in 1990 and is often referred to as Caputo fractional 
derivative: 

1 f* f {v) (o-) 
T{v-a)J a (t-a) 

dt a = V ' 

v - 1 < a < i/, c e N. In contrast to the Riemann- 
Liouville definition it is not necessary to define fractional 
order initial conditions, making Caputo's definition more 
suitable in the context of solving equations with fractional 
derivatives. 

Note that other definitions of fractional derivatives are 
known, like the one due to A.K. Griinwald and A.V. Let- 
nikov which is especially useful when dealing with dis- 
crete approximations (e.g. Schmidt and Gaul (2002)). 
For a wide class of functions the Riemann-Liouville and 
the Griinwald-Letnikov definition are equivalent (see Pod- 
lubny (1999)). 

In literature on parameter identification, most researchers 
seem to use the definition (1) due to Riemann and Liouville 
(e.g. Torvik and Bagley (1984); Bagley (1983); Oldham 
and Spanier (1974); Malti et al. (2007)). However, in most 
of the cases homogeneous initial conditions are assumed 
where (1) and (3) are equivalent. Here, both definitions 
will be treated for homogeneous and inhomogeneous initial 
conditions. 

2.2 Operational calculus 

In this paper Mikusihki's operational calculus is used (see 
Mikusinski (1983)). A brief introduction of this calculus 
in the context of fractional systems can be found in 
Hotzel and Fliess (1998). However, readers unfamiliar with 
this calculus can (in a simplified manner) consider the 
operational expressions as Laplace transforms. 

Using Mikusinski's operational calculus (1) reads 



D?f = 



s a f-Y i ^- a ms v - 1 - 



k=0 



(4a) 



and (3) yields 



C D?f = 



v-l 



s Q /-^/( fe )(0)s a - 1 - fe 



k=0 



(4b) 



v - I < a < v, v £ N, where J^ a f(0) is the limit of 
J"~ a f(t) for t — > 0. Note that both correspondences are 
equivalent in the case of homogeneous initial conditions. 
Hence, for (2) it follows 



Jtf 



s~ a f 



a G 



(5) 



Two fundamental properties of this operational calculus 
are 



[/i/ 2 ](t) = (/i*/ 2 )(t) 



h{o-)h{t-o-)do- (6) 



and 



[/'](*) = -(*/)(*) = -*/(*) 



(7) 



Zero is included in the set R+. 



where t : t i-> t is the identity map and f = j^f denotes 
the derivative w.r.t. s. 

3. THE CASE OF HOMOGENEOUS INITIAL 
CONDITIONS 

Based on the fundamental notions above, in this section, 
the identification problem is addressed for fractional sys- 
tems assuming homogeneous initial conditions. This way, 
all basic ideas can later quite easily be adapted to the 
general case. 

3.1 An introductory example 

A commonly used example is the empirical model 

a(t) = Eve{t) + E 1 D«e{t) (8) 
of a viscoelastic material 3 , also referred to as three param- 
eter generalized Voigt model (e.g. Bagley (1983); Podlubny 
(1999)). The stress a is expressed as a superposition of an 
elastic part E e and a viscoelastic part E\Dfe using a 
fractional derivative of order a e (0, 1) of the strain e. 

The aim here is to identify the parameters E a , Ei, as 
well as the fractional order a from known signals a 
and e. Therefore, using definition (4a) for an operational 
notation, the fractional model (8) is written as 

& = E i + E lS a £, (9) 
where for simplicity a homogeneous initial condition 
e(0) = is assumed. 

In order to obtain an equation that can easily be inter- 
preted and that does not involve any fractional derivatives, 
^ is applied to (9): 

a' = E i' + E 1 (i' + as- 1 e)s a . (10) 
Then a combination of (9) and (10) yields an expression 
without s": 

(i' + as" 1 e)(CT - E i) = i(a' - E Q i'). (11) 
Reinterpreting the expressions as functions of time gives 

(s-kia — is*a)(t) = —a (e*a)(r) dr + aE Q (£*e)(r)dr 
Jo Jo 

(12) 

(cf. (6) and (7)). Note that neither fractional derivatives 
nor derivatives of integer order appear in (12). In order 

3 A material is considered elastic for a = and viscous for a = 1. 



to calculate -Bo and from trajectories t h-> e(t) and 
< i— > cr(t) at least one more (independent) equation in the 
parameters is required. 

Eq. (11) together with the equation obtained by multiply- 
ing it with s~ x yields a linear system for a and aE®: 



"/(*■ 



(e*cr)(r) dr 
-t)(£*ct)(t) cJt 



(e*£)(r) cir 
(i— r)(e*e)(r) rfr 

(e-kicr — t£*o")(r) e?r 



o 
o<.Eq 



(13) 



where double integrals have been replaced by simple ones 
using the general relation 



\fe-i 



y(r fe _i) dT fc _i = J — — — y(r) dr 

(14) 

for fc-times integrals, fe > 1. Eq. (13) can be solved for the 
parameters as long as the matrix in the relation is regular. 
Obviously the matrix is singular for zero signals e and a. 
For all other signals the determinant is expected to vanish 
only at singular points. 

Based on the parameter estimates obtained from (13), 
the remaining coefficient E\ can easily be calculated by 
means of (8). However, since differentiations of e are to be 
avoided, (8) is integrated (corresponding to a multiplica- 
tion of (9) with s _1 ). Solving for the missing parameter 
then gives 

g o{t) dr - E g £ (r) dr 

El ~ ■ (15) 

For validation purposes a simulation of the introductory 
example (8) was performed with values for Polybutadiene 
as in Bagley (1983). The parameters Eq and E\ as well as 
the fractional order a are identified using (13) and (15). 
The integrals in (15) are approximated using 

a N-l 

J2A k+1 f(t-k±) (16) 

fc=0 



Jtf{t) 



with N G N and 
A k +i = 



T(k + a) 



k - 1 



r(a)r(fc + i) k ~ Ak > (17) 

which follows from the Griinwald-Letnikov definition of 
fractional integral and derivatives (cf. Schmidt and Gaul 
(2002)). Fig. 3 shows exemplary results of the identifica- 
tion when using N = 4001 and signals with white noise as 
depicted in fig. 3. For t > 3 s, all parameters are identified 
with an error of less than two percent. The results can be 
improved even further by integrating (13) to reduce the 
impact of the white noise, and by considering additional 
equations and solving a least squares problem. 

3.2 The general approach 

In the case of homogeneous initial conditions the defini- 
tions (1) and (3) of fractional derivatives are equivalent. 
Therefore, without loss of generality the one due to Rie- 
mann and Liouville is used in this section. 



Consider the following relation between two known sig- 
nals 4 u and y: 

(So + ai^ t ai +•••+«„ D?")y(t) 

= (b + hDf" +1 +■■■ + b m D?" +m )u(t). (18) 

At least some of the parameters cij G R, i = 0,...,n 
and hi G K, i = 0, . . . , m as well as the fractional orders 
on G R+, i = l,...,n + m are considered unknown 
and are to be identified. Note that if all coefficients in 
(18) are unknown, only the coefficients in a normalized 
equation can be identified. Define by O the set of unknown 
parameters. 

With (4a), the operational representation of (18) reads 

(a + ais" 1 H h a„s Q ") y 

= (b Q + b 1 s a "+ 1 +--- + b m s a "+™)u. (19) 

First, the fractional operators s Qi , i = 1, . . . , m + n are 
eliminated. For that, (19) is rewritten as 

(a (s)y - b (s)u) S a « + ( ai (s)y - h(s)u)s ai 

+ ■■■ + (a r (s)y~b r {s)u)s a '' =0 (20) 

where ai(s),bi(s) G K(s), i — 0,...,r with r < n + m 
are rational in s. These expressions obviously involve the 
original parameters on, i — 0, . . . , n and bi, i = 0, . . . , m. 



4 The generalization towards multiple equations with potentially 
more than two signals is obvious. 
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Fig. 1. Identified values for a, E , and Ei (solid) compared 
to model parameters (dashed). 
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Fig. 2. Trajectories of stress a and strain e. 



Note that one or more of the fractional orders on, i = 
0, . . . , r can be integer or even zero. 

As in the introductory example, (20) is derived w.r.t. s a 
total of r times. The resulting homogeneous linear system 
of equations for s ai , i = 0, . . . , r has the general structure 

/ a (s)y - b a (s)u ■■■ a r (s)y - b r (s)u \ 

Pl,o(s,s)y-<Zl,o(s,^)" ••• Pl,r(s,^)y-gi, r (s,^)M 

\Pr,ds,j- s )y-qr,ds,-f s )u ••• p r , r (s,-f s )y-q r . r (s,-^)uj 

• ; =o (21) 
W 

where pi j(s,£),qij{s,f s ) G R(s)[£], i = l,...,r, j = 
0, ...,r are polynomials in with coefficients in R(s). 
They are defined recursively by 

p*A s > £) = foi-iA*> j- s )+p^M £)(£ + 

(22a) 

nA*, is) = &Qi-iA*> is) + Qi-iA*> tsKts + a o s ~ 1 ) 

(22b) 

wiihp 0J (s,£) = aj(s) andq 0ij (s,£) = bj(s),j = 0,...,r. 
Note that products p i _j(s,j s )y and q i j(s,^)u can also 
be written as sums of y and u and its derivatives with 
coefficients in R(s). Hence, the operator ^ does not 
explicitly occur. 

From (21) it becomes apparent that the matrix in this 
equation, henceforth denoted by P, has to be singular 5 . 
Then, from det P = an equation is obtained that involves 

• (most importantly) only integer orders of s, 

• products of r operators u, y, and their derivatives 
w.r.t. s, 

• as well as products of parameters a» € R, i = 0, . . . , n 
and hi e R, i = 0, . . . , m and fractional orders 
on e K+, i = 0,...,r. 

A multiplication with s~ k with sufficiently large k such 
that only non-positive powers of s remain gives an expres- 
sion of the form 

F(s-\ y,...,y^\ «,..., «W)=0, (23) 

where is polynomial w.r.t. to all its arguments and 
parameters 0i C 6 including all fractional orders on, 
i = 0, . . . , r to be identified. The knowledge of at least 
one of the orders is required in order to identify the others 
(cf. (20)) - again, this can be considered as a kind of non- 
restrictive normalization. 

Note that the expressions occurring in (23) can easily be 
written as functions of time t again. Hence, depending on 
the number of unknown (independent) coefficients in (23), 
which is larger or equal to the cardinality of 0i, additional 
equations are required to solve for the parameters. They 
are obtained, e.g., by multiplication of (23) with negative 
powers of s. Further integration of the corresponding time 
functions yields independent equations for sufficiently rich 
signals. While other approaches of generating independent 
equations are possible - basically by applying other op- 

5 The matrix P is the result of applying the operators, the entries 
of which belong to a commutative ring. 



erators from R(s)[^] to (23) - which might exhibit a 
better numerical behavior or might be more suitable in 
the presence of noise, this is not in the focus of the present 
contribution. Also, the issue of whether the parameter 
problem arising from (23) is solved as a nonlinear equation 
or a linear one (by overparamctrization) is not discussed 
here. Several approaches are well-known to the community. 
In the context of algebraic identification some insight into 
this is given in Fliess and Sira- Ramirez (2003) and Gehring 
et al. (2012). 

It is rather obvious that parameters 62 C with 6 = 
81 U 02 exclusively occurring as a factor of a fractionally 
differentiated signal in (18) or more precisely a common 
factor of any a,i(s)y — bi(s)u, i = 0, ...,r in (20) will 
not occur in (23). However, as done for the introductory 
example, these parameters can be calculated based on the 
identified values for 9i and (20). 

To this end, based on estimates for the fractional orders, 
(20) is multiplied with s~", N 3 v > max^c^, i.e., such 
that only non-positive powers of s remain. Then, as before, 
depending on the number of parameters remaining to be 
identified (which corresponds to the cardinality of Q2) 
additional equations are obtained by multiplication with 
negative powers of s. The solution of this linear system of 
equations completes the identification. 

4. THE CASE OF INHOMOGENEOUS INITIAL 
CONDITIONS 

For inhomogeneous initial conditions the definitions (1) 
and (3) are essentially different. To illustrate this, first, 
the introductory example is revisited. 

4-1 Introductory example 

Riemann-Liouville fractional derivative For inhomoge- 
neous initial conditions the operational representation of 
(8) (with the Riemann-Liouville operator) reads 

a = E i + E x (s a e - J^ a e(0)) (24) 
since a € (0, 1) and, therefore, v = 1. Again, (10) 
results from applying to (24) and is rewritten here for 
completeness: 

&' = S e' + ^i(e' + as- 1 e)s a . (25) 

If the initial value J t 1_Q e(0) is known, both equations can 
be combined such that s a is eliminated. Similarly as in the 
case of homogeneous initial conditions, (25) is obtained 
with a + EiJ]~ a £(0) instead of a: 

(£' + as~ 1 e)(6 + E 1 j}-' a e{0)-E s) =i(a'-E i'). (26) 

However, since the initial value J t 1_Q e(0) is not simply the 
value of £ at t = but is more complicated to interpret, it 
is only plausible to consider it as unknown. As discussed in 
the sequel, the unknown initial values (here only one) can 
either be eliminated by some algebraic manipulations or 
they can be treated as additional parameters and included 
in the identification. For the present example, the latter 
approach simply means that there are three parameters 6 - 
or products of parameters - namely a, E a , and E\ J t _a £(0) 
to be identified in (26). 

6 The remaining parameter Ei is again calculated from (15). 



On the other hand, an elimination of an unknown initial 
value bears the advantage of being faced with less param- 
eters to be identified. Clearly, (25) does not depend on 
e(0) anymore. So deriving once more w.r.t. s allows for 
the elimination of s a : 

(£" + 2as- 1 e' + a(a - l)s~ 2 £) (a' - E i') 

= (£' + as- 1 ^" - E £"). (27) 
The operational expressions in this equation can again 
be written as functions of time, which yields a nonlinear 
equation in E and a. In order to determine the unknown 
parameters either the nonlinear parameter problem can 
be solved or overparametrization can be used treating the 
coefficients a, a 2 , E , aE , and a 2 E as independent. 

Caputo fractional derivative The approach is a little 
different if Caputo's fractional derivative is used instead of 
the one due to Riemann and Liouville. In order to illustrate 
this, Df in (8) is replaced by C D" . The correspondence 
(4b) then yields 

a = E o e + E 1 (s a e-e(0)s a - 1 ). (28) 

As mentioned before, Caputo's definition bears the advan- 
tage of using initial values (here e(0)) that can easily be 
interpreted. While in this simple example the presumed 
knowledge of e(0) does not impose any restriction, since 
e itself is considered as a known signal, for a > 1 initial 
values of derivatives come into play that are unknown in 
general. For that reason, and because in a real application 
measurements arc noisy, e(0) is assumed to be unknown 
in the sequel. 

Deriving (28) w.r.t. s gives 

a' = E e , +E 1 (s a (e'+as- 1 e)-(a-l)s- 1 e(0)s a - 1 ). (29) 

If the initial value is to be identified, eliminating s a from 
(29) using (28), 

{&'£ - &£') = (e(0)(l - a)s~ 2 + as^e)^ - E e) 

+ e(0)s- 1 (cr' - Eoi') (30) 

is obtained which is an equation of the type (23) and 
involves the parameters e(0), E 0} and a. 

Alternatively, the initial value e(0) can simply be elimi- 
nated by multiplying (28) with (1 — a)s _1 and adding it 
to (29) multiplied with s: 

(l-a)a + s&' = E (l-a)i + E s£' + E 1 s a (i + si'). (31) 

By taking yet another derivative of this equation w.r.t. s 
the factor s a corresponding to the fractional derivative is 
replaced. Multiplying with s~ 2 , the resulting equation 

a'e" - a"£' + S - l {&£" - a"e) + 2 S - 2 {ae' - &'£) 

= as- 2 (ae' - &'£) - as- 1 {2£' + s^e)^' - E £') 

+ as- 1 (£" + as- 1 ? + (a - l)s~ 2 e) (a - E e) (32) 

is of the type (23) and involves the parameters E and a. 

4-2 The general approach 

Based on the previous example and the discussions in 
section 3.2 for the homogeneous case the general approach 
for an equation (18) is rather obvious. That is why only a 
sketch is given of how to obtain relations of the form (23) . 
For both the Riemann-Liouvillc and the Caputo definition 
the cases where initial conditions are eliminated and the 



one where they are identified are treated. Henceforth, it is 
assumed that an upper bound of all unknown (integer or 
fractional) orders on, i = 1, . . . , n + m is available 7 . 

Riemann- Liouville fractional derivative For a known 
upper bound v e N of the order of differentiation, the 
operational representation of (18) (corresponding to (20) 
for homogeneous initial conditions) can be written as 

(a (s)y - b (s)u)s a ° + ( ai (s)y - h{s)u)s ai 

+ ■■■ + (a r (s)y - b r (s)u)s^ = £- =1 as*' 1 (33) 

where the coefficients Cj, i — 0,...,v depend on initial 
values of u and y and their derivatives as well as parame- 
ters ai, 0, . . . , n and bi, i = 0, . . . , m. Applying ^ to this 
equations eliminates all these initial values. Following the 
elimination of the fractional operators s ai , i = 0, . . . , r, a 
multiplication with s~ fc (with sufficiently large k) yields 
an equation of the form (23) with r + v instead of r. 

On the other hand, identifying the initial values - even 
though they lack physical meaning - is almost identical 
to the procedure in section 3.2, the only difference being 
that (21) is inhomogeneous with the right hand side of 
the equation being a vector of polynomials in s. Using 
the solution of this linear problem and substituting it 
in an additional equation involving the fractional orders 
s ai (generated from deriving once more w.r.t. s), again, a 
multiplication with a sufficiently large negative power of s 
results in (23). 

Caputo fractional derivative The operational represen- 
tation of (18) with Caputo's derivation operator gives 

(a (s)y - b {s)u - Yh=i c Q . t s- l )s aa 
+ ■■■ + (a r (s)y - b r (s)u - Cr,i* _i )* ap = ( 34 ) 

for a known upper bound v € N. The coefficients Ci_j, 
i = 0, . . . , r, j — 1, . . . , v can be cancelled by taking a 
sufficient number of derivatives w.r.t. s and solving the 
resulting linear system of equations for these coefficients. 
Afterwards, an expression of the form (23) can be obtained 
using r further differentiations of this result and multiply- 
ing it with a sufficiently large negative power of s. 

If, however, one is interested in the initial values, as in 
the homogeneous case, the fractional operators s ai in (34) 
can simply be eliminated by taking r derivatives and 
calculating the solution of this linear system (cf. (21)). 

5. LINEAR PARTIAL DIFFERENTIAL EQUATIONS 
WITH FRACTIONAL DERIVATIVES 

The identification of fractional orders and model param- 
eters discussed for lumped fractional models can be ex- 
tended to fractional distributed parameter systems. For 
the integer order case, Rudolph and Woittennek (2008) 
and Gehring et al. (2012) demonstrated that based on the 
knowledge of the solution of the corresponding operational 
ordinary differential equation a relation of the form (23) 
between two boundary measurements can be obtained. 

Here, the basic ideas of the identification method are 
demonstrated on an example. Since homogeneous initial 

7 This assumption constitutes only a mild restriction since any 
(sufficiently large) finite number could be used as the upper bound. 



conditions are assumed, without loss of generality, only 
the Ricmann-Liouvillc fractional derivative is used. 

The example considered is the fractional diffusion-wave 
equation (see e.g. Mainardi (1997); Podlubny (1999)) 

v 2 D 2 z u(z,t) = D?u(z,t), (z,t)eR 2 + (35) 

with Dgl, where the derivative w.r.t. time is of fractional 
order a G (0,2]. The factor v 2 is the diffusion coefficient 
for a = 1, and it is the squared wave propagation speed for 
a = 2. In the context of signal propagation as addressed 
in Mainardi (1997) the boundary conditions are 

u(0,t) = h(t), u(oo,t) = 0. (36) 

Additionally, as mentioned before, homogeneous initial 
conditions are assumed 8 . 

The operational representation of (35) reads 

v 2 D 2 z u{z) = s a u(z) (37) 

and is a second order ordinary differential equation w.r.t. 
z that can easily be solved. From the general solution 

u(z) = C^""'* + C 2 e-^ sa/2 (38) 
using the (operational form) of the boundary conditions 
(36) yields operators C\ = and C 2 = h. 

In order to identify the fractional order a, apart from 
the knowledge of h, another measurement is required. 
If u is measured at some point z = L =/= 0, such that 
g(t) = u(L, t) is known, (38) gives 

g = he-^ sa/2 . (39) 

Since the exponential function satisfies a first order differ- 
ential equation it can be eliminated by deriving (39), 

g' = {h! - 1 ^s^s^h) e"^ s ° /2 (40) 

and replacing the exponential function using (39): 

gh' -g'h=%±s- 1 s a / 2 gh. (41) 

Finally, s a / 2 is eliminated using the derivative of (41) 
w.r.t. s. The interpretation of 

g 2 ((h'f s-'hh' hh") h 2 ((g'f s-'gg'-gg") 

= ^8- 1 ghtfh-gh') (42) 

in terms of functions of time is omitted for brevity. It in- 
volves nested convolutions of known signals only. However, 
as can clearly be seen, (42) constitutes a linear equation 
in a. Based on an estimated value for the fractional order, 
the quotient 9 L/v can directly be calculated from (41) 
(since s a / 2_1 already corresponds to a fractional integral 
for a e (0,2)). 

6. CONCLUSION 

Using the algebraic framework presented model param- 
eters and fractional (as well as integer) orders can be 
identified in linear fractional models. The exact relations 

8 Note that the number of initial conditions depends on the (un- 
known) fractional order a. However, this fact is of no interest in the 
present discussion. 

9 Prom a physical perspective it seems plausible that either the 
distance L or the "wave propagation speed" v have to be known 
in order to calculate the other one (based solely on the knowledge of 
u(z,t) at two points z). 



obtained for the unknown quantities involve only convolu- 
tions of known (input and output) signals. Based on pos- 
itive results using simulation data, the authors intend to 
further validate the identification algorithms using exper- 
imental data. In this context numerical aspects related to 
online implementation and robustness w.r.t. measurement 
noise should be investigated. 

The example of the fractional diffusion-wave equation 
demonstrates, that the method is also applicable to (at 
least certain) linear fractional distributed parameter sys- 
tems as well as delay systems with fractional derivatives. 
Beyond that, an extension towards the identification of 
the structure of linear systems, i.e., ordinary differential 
equations, seems possible where, both, the order of the 
equations and their coefficients can be estimated. 
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